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The charge susceptibility of the two-dimensional repulsive Hubbard model is investigated using 
the diagram technique developed for the case of strong correlations. In this technique, a power 
series in the hopping constant is used. It is shown that once the Fermi level crosses one of the 
Hubbard subbands a sharp peak appears in the momentum dependence of the static susceptibility. 
With further departure from half- filling the peak transforms to a ridge around the F point. In 
the considered range < 1 — n\ < 0.2 of the electron filling n the static susceptibility is finite 
which points to the absence of the long-range charge ordering. However, for 1 — n| ~ 0.12 the 
susceptibility maxima are located halfway between the center and the boundaries of the Brillouin 
zone. In this case an interaction of carriers with tetragonal distortions can stabilize the charge 
density wave with the wavelength of four lattice spacings, as observed experimentally in the low- 
temperature tetragonal phase of lanthanum cuprates. In the range of parameters inherent in cuprate 
perovskites the character of the susceptibility evolution with n depends only weakly on the ratio of 
the nearest-neighbor hopping constant to the Hubbard repulsion and on details of the initial band 
structure. The location of the susceptibility maxima in the Brillouin zone is mainly determined by 
the value of n. 
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I. INTRODUCTION 

A variety of numerical methods have been used to elu- 
cidate the existence of phase separation in the ground 
states of the Hubbard and the related t-J models. The 
states with charge density waves (CDW, stripes) were 
obtained in Refs. 1,2.3.4.5.6 using the mean-field approx- 
imation, the variational principle with the Gutzwiller- 
type variational functions and the density matrix renor- 
malization group calculations. However, results of Monte 
Carlo simulations 7 *^ and cluster calculations! *!! have 
cast doubt on this finding. Hence different techniques 
give different results which reflects the fact that nearly 
degenerate low-lying states of the models have different 
nature. 

Experimentally static stripes were observe d 12 ' 13 ' 14 ' 15 
in the low-temperature tetragonal (LTT) phase 
of lanthanum cuprates Lai.6-^Ndo.4Sr x Cu04 and 
La2_a;Ba x Cu04. One of the manifestations of the 
stripe formation is the anomalous suppression of su- 
perconductivity near the hole concentration x — i in 
La2_ a; Ba 2 ;Cu04. A weaker suppression of T c near this 
hole concentration is also observed in I^-^Sr^CuC^ in 
the low-temperature orthorhombic (LTO) phase.— The 
mentioned phases are characterized by tilts of the CuOg 
octahedra about axes (LTT) and diagonals (LTO) of the 
Cu-0 planes. These experimental observations suggest 
that the interaction of carriers with respective phonons 
plays an essential role in the stripe stabilization. 

For investigating instabilities of a system it is conve- 
nient to use corresponding static susceptibilities. Similar 
to the divergence of the magnetic susceptibility, which 



points to the establishment of the long-range magnetic 
order while finite maxima indicate short-range ordering, 
a divergence of the charge susceptibility means the ap- 
pearance of the CDW, while finite peaks are manifesta- 
tions of respective charge fluctuations. In this paper we 
calculate the charge susceptibility of the two-dimensional 
repulsive Hubbard model using the strong-coupling dia- 
gram technique. In this technique ; 17 ' 18 ! 19 ' 20 on the as- 
sumption of strong electron correlations Green's func- 
tions are calculated using the expansion in powers of the 
hopping constant. The terms of this expansion are ex- 
pressed by means of site cumulants of electron creation 
and annihilation operators. 

It is found that once the Fermi level crosses one of the 
Hubbard subbands a sharp maximum appears in the mo- 
mentum dependence of the static susceptibility near the 
T point. With further departure of the electron filling 
n from the half-filling value h — 1 the maximum trans- 
forms to a ridge around the T point. In the considered 
range of the electron filling < |1 — n\ < 0.2 the static 
susceptibility is finite which implies the absence of the 
long-range charge ordering. However, obtained results 
give some insight into the way in which phonons can sta- 
bilize stripes. For 1— n| ps 0.12 the susceptibility maxima 
are located halfway between the center and the bound- 
aries of the Brillouin zone. In this case an interaction of 
carriers with the tetragonal distortions gives a maximal 
energy gain for the CDW with the wavelength of four 
lattice spacings, as observed experimentally in the LTT 
phase of lanthanum cuprates. The calculations were car- 
ried out for the t-t'-U Hubbard model for the ratios of the 
Hubbard repulsion to the nearest-neighbor hopping con- 
stant U/\t\ = 8 and 12, and for the next-nearest-neighbor 
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hopping constant t' = and — 0.3i. The above-discussed 
evolution of the susceptibility with n is not changed qual- 
itatively with the variation oiU/t and t'/t. The shape of 
the susceptibility and the location of its maxima in the 
Brillouin zone are mainly determined by the value of h. 
The mentioned values of the parameters belong to the 
parameter range which is widely believed to be suitable 
for cuprate perovskites. 

Main formulas used in the calculations are given in the 
following section. The discussion of the obtained results 
and their comparison with Monte Carlo simulations are 
carried out in Sec. III. Concluding remarks are presented 
in Sec. IV. 

II. MAIN FORMULAS 

The Hubbard model is described by the Hamiltonian 
H = ^ tgiC^aVa + Tj- n\ a n\- a , (1) 

ll'cr la 

where a\ a and ai CT are the electron creation and annihi- 
lation operators, 1 labels sites of the square plane lattice, 
a =f or I is the spin projection, t\y and U are hopping 
and on-site repulsion constants, and n\ a — a\ a a\ a . Be- 
low we consider the case where only the nearest-neighbor 
t and next-nearest-neighbor t' hopping constants are 
nonzero. 

As mentioned above, in the considered case of strong 
correlations, U \t\, t' , we use the strong-coupling 
diagram techniqu e 17 ! 18 ' 19 ! 20 for calculating the charge 
Green's function 

B,v(lV,lr) = (TSn Va ,(T')Sn la (r)), (2) 

where Sn\a- — n\ a — (n\ a ) is the deviation of the elec- 
tron occupation number from its mean value, the angular 
brackets denote the statistical averaging with the Hamil- 
tonian 

H = H — jU^JniCT, 

lcr 

\x is the chemical potential, T is the time-ordering oper- 
ator which arranges other operators from right to left in 
ascending order of times r. The time evolution of oper- 
ators in f2} is also determined by the Hamiltonian TC, 

0{t) = exp(Hr)Oexp(-?ir). 

Using the strong-coupling diagram technique one can 
convince oneself that equations for B are similar to those 
derived for the spin Green's function in Ref. [2l|. After 
the Fourier transformation these equations read 

B a , a (q) = -S^Y.Gfa-Wl + Pi) 
Pi 

+ (1?) £ n (Pi) I I(pa)n(«+Pi)n(g + i> 2 ) 

X K><r(pi,q + Pl,q + P2,P2), (3) 



= K'<r{pi,q + pi,q + P2,P2) 

-jr z2 K'at(pi,q + Pi,q+ P3,P3)Q(p3)Q(q + P3) 

xA CTlCT (p3,g + p 3 ,g + p 2 ,P2). (4) 

Here the combined indices q = (k, iu) v ) and pj — 
(k^,io; n .) were introduced, lu u = 2vttT and uj n = (2n + 
l)nT are the boson and fcrmion Matsubara frequencies 
with the temperature T, k is the wave vector, G{p) = 
{( a ko-\ a U a }) is the electron Green's function, H(p) = 
1 + tkG(p), tk is the Fourier transform of the hopping 
constants which is equal to = 2t[cos(k x ) + cos(k y )] + 
4i' cos(k x ) cos(fc y ) in the considered case, the lattice spac- 
ing is taken as the unit of length, 0(p) = t\JI(p) is the 
renormalized hopping, A. l7 i (7 (pi,p + pi,p + P2,P2) is the 
sum of all four-leg diagrams, \ a > a{pi, p + pi, p + P2, P2) is 
its irreducible subset, and N is the number of sites. 

The main difference between Eqs. (|3]), (j4} and the 
respective equations for the spin Green's functional is 
in the spin indices of the irreducible four-leg diagrams. 
For the transversal spin Green's function considered in 
Ref. l2ll only \n(pi , p+Pi, P+P2, P2) enters into the equa- 
tions. 

Equations ([3]) and (|4]) can be somewhat simplified if we 
take into account the invariance of Hamiltonian |T]) with 
respect to rotations of the spin quantization axisi 2 ^ This 
invariance leads to the following symmetry relations: 

B n (q)=B u (q), B u {q) = B n (q) 

and analogously for A CT v and \ a ' a - Using these relations 
we find from Eqs. ([3]) and Q 

pi 

+ (jf) En( Pl )n(p 2 )n(g + Pl )n( (? + p 2 ) 

^ ' P1P2 

X A(p 1 ,q + p 1 ,q + p 2 ,P2), (5) 
k(pi,q + Pi,q+P2,P2) = HPi,q + Pi,q + P2,P2) 

_ jy H 1 + P 1 ' 1 + P3,P3)@(P3)9(g + Pi) 

P3 

xA(p 3 ,q + p 3 ,q + p 2 ,P2), (6) 

where 

a' (7 

A(p 1 ,q + p 1 ,q + p 2 ,p 2 ) 

= 2 ^2K'a(pi,q + Pi,q + P2,P2) (7) 

a' a 

X(p 1 ,q + p 1 ,q + p 2 ,P2) 

= 1^^2K'a(Pl,q+Pl,q+P2,P2)- 
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FIG. 1: Diagrams corresponding to Eqs. ([5]) and (B). 

Diagrams corresponding to Eqs. (O and ([6]) are plot- 
ted in Fig. Q] In these diagrams, the dual line indicates 
the electron Green's function G(p), the shaded and open 
squares stand for A and A, respectively, the single di- 
rected line between the squares is the renormalized hop- 
ping 0(p), and the dashed line is the external line of the 
four-leg diagram with the inserted irreducible two-leg di- 
agrams. The dashed line corresponds to II(p). 

In the following consideration, we simplify the general 
equations (0 and © by neglecting the irreducible two- 
leg diagrams in the external and internal lines of the four- 
leg diagrams. In this approximation quantities II(p) and 
8(p) in Eqs. ([5]) and ([6]) are substituted by 1 and ik, 
respectively. Besides, we use the lowest-order irreducible 



four-leg diagram instead of \r' a {pi,q + pi,q + p 2 ,p 2 )- 
This four-leg diagram is described by the second-order 
cumulant of electron operators 

K%'°(T',T,T[,n) = (Ta CT ,(r> CT ,(T)a CT KK(n)}o 
-Ki(t', t)K\ (t[,ti) 

+K l {T l ,T l )K l {T l ll T)5 a , a , (8) 

where the subscript of the angular brackets indicates 
that the averaging and time dependencies of operators 
are determined by the site Hamiltonian 

Hi = Y^[(U/2)ni <7 ni,- (T - pn Xa ], (9) 

u 

Sio-(t) = exp(H\T)al a exp(— Hit), and the first-order cu- 
mulant Ki(t',t) = (Ta (T (r')a cr (T)}o. This cumulant 
does not depend on the spin index a and the respec- 
tive superscript is dropped in the cumulant notation. All 
operators in the cumulants belong to the same lattice 
site. Due to the translational symmetry of the problem 
the cumulants do not depend on the site index which 
is therefore omitted in the above equations. After the 
Fourier transformation the expression for K 2 a reads 



K. 



2 (ni,v + m,v + n 2 ,n 2 ) 



^{^[(XoiW -S nin2 )e- El13 
+Z~\8 vQ <W 2 <W) (e-^°+^ - e- 2E ^)] F( ni + v)F{n 2 ) 

-8 a >^ a e~ Bol3 UgQi(ni + v)gm{n 2 )gn 2 {n 1 + n 2 + v)[goi(n 2 +v) + .9oi(«i)] 
-Sa'-ue'^PUgx^nx + v)gi 2 {n 2 )g Q2 {n 1 + n 2 + v) [gi 2 {n 2 + v) + .gi 2 (ni)] 



+Sa'. 



F(m + v)goi(n 2 )g i(n 2 + v) + F(n 2 )g m (ni + v)g m (ni 



-F(n 2 )g 12 (n 2 + v) [gx 2 {ni + v) -goi(ni)] + F{nx + v)g\ 2 {ni) [gi 2 {n 2 ) - g i(n 2 + v)] j, 



(10) 



where Eq = 0, E\ = —p, and E 2 = U — 2/j, are the 
eigenenergies of the site Hamiltonian (3 — T^ 1 , Z — 
e -E f3 _|_ 2e~ Elf3 + e~ E2t3 is the site partition function, 
9ij( n ) = H + Ei - Ej)- 1 , F(n) = g i(n) - gu(n), 
and integers n and v stand for the fermion and boson 
Matsubara frequencies. 

Equation (|10p can be significantly simplified for the 
case of principal interest U ^> T. In this case, if pL satisfies 
the condition 



r) < n <U -77, 



(11) 



where ij ^> T, the exponent e @ El is much larger than 
e -/3£o an d e -P E 2. Therefore terms with e _/3Bo and e~ pE2 
can be omitted in Eq. (|10|) . In this case we obtain for the 
quantity K 2 — i ^2 a i a K% a which is used instead of A in 



Eq. © 

K 2 (ni, v + m, v + n 2 ,n 2 ) 
3 

= --l38 nin2 F{m + v)F{n 2 ) 

+ -^F(n 1 + v)goi{n 2 )g Q1 (n 2 + v) 

+F(n 2 )g 01 (n 1 + v)g al {n 1 ) 
+F(n 2 )g 12 (n 2 + v) [p l2 (ni + v) - g 01 (n x )] 

+F{n\ + v)gi2{n\)[g\ 2 {n 2 ) - goi(n 2 + v)] }■ (12) 

As can be seen from Eq. (|6|), in this approxima- 
tion A(k x , u ni ; k + ki , oj u +uj ni ; k + k 2 , w„ +u n2 ; k 2 , w„ 2 ) 
does not depend on ki and k 2 . 

Taking into account these simplifications we find for 



4 



the second term on the right-hand side of Eq. ([5]) 
j 1 \ 2 g 

— J ^ A = --T^/ k (im)ai(2/ + n)ai(ri) 
' P1P2 n 

~2 [1 + S'kjMki')]! 1 ^ fk(vn)a 3 (-i/, v + n) 

n 

- S'kjyi (ki/)] r ^ h{vn)a 2 (—v, v + n) 



ryj/k(tm)o4(-i / i v + n) 

n 

7;S k y 2 (k.v)Ty^f k (L>n)ai(i' + n), 



(13) 



where 5 k = N 1 X) k , *k+k'*k', 



1 - -S'ki ? (^ + n) J F 1 (n) 



/k(fn) 

ai(n) = F(n), a 2 (vn) = g 01 (n)g 01 (v + n), 



a 3 (vn) = F(n)g 12 (v + n), 
a 4 (im) = 3i 2 (n) - 501 + ?U, 



(14) 



and yiQsp) are solutions of the following system of four 
linear equations: 

2/i(kf) 

= T 2 X! /k( (y "-i)^2("i, ^ + ni, v + n 2 , n 2 )a i (vn 1 ) 



\^ kT ^l /k(fn)ai(i/ + n)cij (1/71)3/2 (ki/) 
^S k T^/ k (i/n)a 2 (-^ z/ + n)ai{vn)yi (kz/) 

n 
n 

-^ s kT~^2 fk(m)a 3 (-v, v + n)ai(vn)y 4 (kv). (15) 



The system of equations (|T5J) follows from Eq. ([6]), 
Vi(kv) — T 2 2J ai(^i)Ak(«i, v + f + n 2 , n 2 ). 



In calculating the first term on the right-hand side 
of Eq. ([5]) we use the Hubbard-I approximation^ for 
the electron Green's function. Under condition (fTTj) this 
function reads 



G(kn) 



iu; n + fi — U/2 



(iu) n + [i){iuj n + n - U) - t u (iuj n +fi- U/2) ' 

(16) 
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FIG. 2: The static charge susceptibility in a 4x4 lattice for 
t = -U/8, t' = 0, and T = 0.125|t|. The susceptibility is plot- 
ted along the triangular contour in the Brillouin zone. The 
corners of the contour are given by the momenta k = (0, 0) 
(r), (7r,0) {X), and (w, 7r) {M). Open symbols are results of 
Monte Carlo simulations 24 for n = 1 (squares) and n = 0.95 
(circles). Filled symbols are our results for the same electron 
fillings. 



III. CHARGE SUSCEPTIBILITY 

In this section we consider the static charge suscepti- 
bility, 

Xo (k) = B(k,i/ = 0). 

As mentioned above, divergence of this quantity points 
to charge instability of the system, while its finite max- 
ima indicate regions of increased charge response on an 
external field. To check the used approximation we com- 
pare our results obtained for a 4x4 lattice with the data 
of Monte Carlo simulations carried out for the same 
conditions.— This comparison is demonstrated in Fig. [2] 
In this figure, n — ^2 a {ni a ) is the electron filling. The 
value n — 1 corresponds to half-filling. We used a small 
global vertical shift of our results to fit the Monte Carlo 
data better. As seen from the figure, our approximation 
reproduces correctly the general shape of the dependence 
Xc(k) and its variation with electron filling. Notice the 
rapid increase of the uniform susceptibility Xc(0) with 
departure from half-filling. 

The susceptibility for a larger lattice and a smaller 
temperature is shown in Fig. [3J At half-filling the sus- 
ceptibility is small and its dependence on momentum is 
weak, as can be seen in Fig. EJa). This shape of Xc(k) is 
retained, until the Fermi level crosses one of the Hubbard 
subbands. Immediately after this crossing, which leads 
to departure from half-filling, a sharp peak appears in 
the dependence Xc(k). As shown in Fig. [31(b), the peak 
is finite and is located near the T point which indicates 
the appearance of long-wave charge fluctuation in the 
system. With increasing x = |1 — n\ the susceptibility 
grows. For moderate values of x it peaks for momenta 
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along a contour which is centered on the T point and has 
the shape of a somewhat distorted circle [see Fig. [3][c)]. 
With increase of x the size of the contour grows indicat- 
ing the decrease of the wavelength of charge fluctuations. 
Notice that in Fig.[3Jc), for x « 0.12, the contour crosses 
the axes at k Xl k v rj tt/2. On approaching the boundary 
of the considered range of /i, Eq. (TlTj) . the susceptibil- 
ity start to grow near the X and Y [(0, 7r)] points [see 
Fig. EJd)] . For chemical potentials in this range the con- 
tour does not cross the boundary of the magnetic Bril- 
louin zone - the line segment connecting the X and Y 
points. Notice also that in this range corresponding to 
electron fillings < x < 0.2 which are most interesting 
for cuprate perovskites the susceptibility remains finite. 
This means that a long-range or stripe ordering of charges 
does not occur in the Hubbard model for the considered 
parameters. The model produces only charge fluctua- 
tions with wave vectors corresponding to the maxima of 
the susceptibility. The susceptibility varies only slightly 



with a further decrease of temperature from the value 
used in Fig. [3] 

The results shown in this figure were obtained for the t- 
U Hubbard model with the ratio of parameters U/\t\ = 8. 
To check how this overall picture is changed with the 
variation of the parameters we carried out analogous cal- 
culations for the t-U model with U/\t\ = 12 and for the 
t-t'-U model with U/\t\ = 8 and t' = -0.3t. This lat- 
ter ratio of the hopping constants is believed to describe 
adequately the band structure of cuprate perovskites. 25 
In both cases the evolution with filling and the shape 
of the susceptibility look qualitatively the same as those 
shown in Fig. [3] Moreover, the location of the suscep- 
tibility maxima in the Brillouin zone varies only slightly 
with the change of the hopping and repulsion constants 
and is mainly determined by the electron filling. As an 
example the susceptibilities for the two considered sets of 
parameters are shown in Fig. [3] for n « 0.88. These sus- 
ceptibilities should be compared with Fig. [3][c). Notice 
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FIG. 4: (Color online) The static charge susceptibility in a 
40x40 lattice for t = -(7/12, t' = (a), and for t = -U/8, 
t' = -0.3t (b). In both cases n « 0.88 and T = 0.001(7. 

that, as in Fig.^c), the maxima on the axes of the Bril- 
louin zone are located approximately halfway between its 
center and boundary. It is worth noting also that in the 
t-t'-U model, which does not possess the electron-hole 
symmetry, the shapes of the susceptibility are qualita- 
tively similar for fillings n = 1 — x and 1 + x. 

The appearance of maxima in the susceptibility with 
departure from half-filling is connected with fermion 
poles of the function A, specifically with the poles given 
by the function f k (vn) in Eq. HH). In general these poles 
do not coincide with poles of the electron Green's func- 
tion G(p). For v — we find four poles of the function 
j\Sv — 0, z) at frequencies 



£< (k) = -/i + | ± J ± ^ llPS*. (17) 

After the transformation to real frequencies sums over n 
in Eqs. (|13p and (TlS)) are transformed to sums over the 



above poles and the terms of these sums contain multi- 
pliers n F (si(k)) where np(a;) = [exp(/3w) + 1] _1 is the 
Fermi-Dirac distribution function. Thus, for low temper- 
atures and small to moderate departure from half-filling 
these sums and the second term of B(q), Eq. (fT"3]) , will 
have a step-like shape as functions of momentum if one 
of the dispersions (fT7|) crosses the Fermi level. This be- 
havior is demonstrated in Figs. G2b), [21c), and|3J 

The above results do not support a purely electronic 
mechanism of long-range charge ordering. However, they 
give some insight into the way in which phonons can sta- 
bilize stripes. As mentioned above, an essential role of 
certain CuOe octahedra tilts in such stabilization can 
be supposed from experimental result s 12 ' 13 ' 14 ' 15 on lan- 
thanum cuprates. It is know n 26 ' 27 that such tilts are 
strongly coupled to the carriers. This interaction can 
be described using an effective electrostatic potential (j>\ 
which is connected with the charge fluctuations {5n\) 
arising in Cu-0 planes as a consequence of the tilts. 
These fluctuations give the following contribution to the 
adiabatic potential: 

AS = -\ J2{Sni)ct>i = -\ E Xc(k)0t (18) 

1 k 

Here we took into account that the charge susceptibility 
is the response of the system on the external potential. In 
Eq. (p~8|) we allowed also for low frequency of the octahe- 
dra tilts.— The potential (f> k depends on the coordinates 
Qkx, Qky of these vibrations. Due to the symmetry of the 
problem with respect to a reflection in a Cu-0 plane this 
dependence contains only even powers of the coordinates 
and starts from quadratic terms, 

0k = A k Ql x + B k Ql y + C k Q kx Q ky + ... (19) 

Hence the contribution (| 18[) is a quartic function of the 
coordinates. As follows from Eqs. (fl~8|) and (fT9|) . finite 
values of the vibration coordinates give an energy gain 
and the larger the value of the susceptibility is, the larger 
energy gain can be achieved (the lattice stability is pro- 
vided by other quartic and higher power terms of the 
adiabatic potential). For weak momentum dependen- 
cies of other parameters of the adiabatic potential the 
largest energy gain is achieved for momenta at which 
the susceptibility is peaked. From symmetry consider- 
ations it can be supposed that for the LTT distortions 
{Qkx 7^ 0, Qky = and vice versa) this gain is achieved 
for maxima on the axes of the crystal plane, and for the 
LTO distortions (Qkx = ±Qky) on its diagonals. Based 
on the experimental result s 12 ' 13 ' 14 ' 15 it is clear that for x 
in the range 0.12 ± 0.02 and for temperatures T < 50 K 
(the domain of the LTT phased in L^-^BaajCuO^) pa- 
rameters of the adiabatic potential provide the lowest 
minimum for tetragonal distortions, and outside of this 
region for orthorhombic distortions. As can be seen in 
Figs. EJc) and [31 for x « 0.12 the susceptibility max- 
ima are located on the axes approximately halfway be- 
tween the r and X (Y) points. Such maxima give the 
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lowest energy for the CDW with the wavelength equal 
to four lattice spacings, as observed experimentally in 
the LTT phase of lanthanum cuprates . 12 d 3 i 14 i 15 Such a 
CDW is commensurate with the lattice which is essen- 
tial for its stability. On the contrary, for x « 0.12 an 
oblique CDW in the LTO phase is in general incommen- 
surate [see Figs. EJc) and 0] which can explain its more 
disordered character J£ 

Another conclusion which can be made from the anal- 
ysis of experimental dat a 14 > 15 and from the above results 
concerns the magnetic incommensurability observe d 29 ^ 30 
in p-type cuprate perovskites in a wide range of hole con- 
centrations, namely that it is unlikely that the magnetic 
incommensurability is connected with charge stripes. In 
La2- K Ba K Cu04 the LTT phase with stripes along Cu-0 
bonds exists in a narrow range of hole concentrations near 
x = 0.12.— Outside of this range, in LTO phase, stripes 
are more disordered 1 '' and oblique. Based on the assump- 
tion that magnetic incommensurability is connected with 
the stripes it is difficult to explain the location of the 
low-frequency incommensurate maxima in the magnetic 
susceptibility along the axes of the Brillouin zone in both 
phases and a smooth variation of the incommensurabil- 
ity parameter with doping through these phases. Mech- 
anisms of the magnetic incommensurability in strongly 
correlated systems, which are not based on charge stripes, 
were considered in Refs. I2lll3lll32ll33l . 

IV. CONCLUDING REMARKS 

In this paper we calculated the charge susceptibility of 
the two-dimensional repulsive t-t'-U Hubbard model us- 
ing the strong-coupling diagram technique. In this tech- 
nique, on the assumption of strong electron correlations 
Green's functions are calculated using the expansion in 
powers of the hopping constant. The terms of this expan- 
sion are expressed by means of site cumulants of electron 
creation and annihilation operators. For a small lattice 
we found a good agreement of the results obtained in this 



approach with Monte Carlo simulations. In a larger lat- 
tice it was found that once the Fermi level crosses one 
of the Hubbard subbands a sharp peak appears in the 
momentum dependence of the static charge susceptibility 
near the T point. With further departure from half- filling 
the peak transforms to a ridge around this point. These 
results were obtained for three sets of parameters: for the 
t-U model with ratios of the parameters U/\t\ = 8 and 12, 
and for the t-t'-U model with U/\t\ = 8 and t' = -OM. 
We found that the shapes of the susceptibility and the 
positions of its maxima are similar in all three cases and 
are mainly determined by the value of the electron filling 
n. This suggests that the obtained shapes of the suscepti- 
bility are inherent in the parameter range of cuprate per- 
ovskites. In the considered interval of the electron filling 
< |1 — n| < 0.2 the static susceptibility is finite which 
points to the absence of the long-range charge ordering. 
However, obtained results give some insight into the way 
in which phonons can stabilize stripes. For |1— n| ~ 0.12 
the susceptibility maxima are located halfway between 
the center and the boundaries of the Brillouin zone. In 
this case an interaction of carriers with the tetragonal 
distortions gives a maximal energy gain for the charge 
density wave which is oriented along the Cu-0 bond and 
has the wavelength equal to four lattice spacings. Such 
stripes are observed in the low-temperature tetragonal 
phase of lanthanum cuprates. We supposed also that the 
obtained incommensurate location of the susceptibility 
maxima along diagonals of the Brillouin zone can be the 
reason of a more disordered character of oblique stripes 
in the low-temperature orthorhombic phase in compar- 
ison with commensurate stripes of the low-temperature 
tetragonal phase. 
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